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I. PURPOSE 


This note describes the methodology of the ZOT,14 program for the 
simulation of the post-enablement phase of flight of a cannon-launched 
guided projectile (CLGP). The mathematical models employed and their 
implementation in the program are discussed. A user's guide is appended 
(Appendix A), 


II, GENERAL DESCRIPTION 


The guidance simulation model ZOT.14 is a general-purpose, quick- 
running, continuous-time simulation model for the post~enablement phase of 
the flight of a guided projectile with proportional navigation, It 
incorporates a non~rolling-projectile flight model (which accounts in large 
part for its speed) with simplified flight equations, a model for the random 
motion of the designator spot about the target, a model for the evasive 
maneuver of the target, and a model for acquisitton/loss of acquisition of 
the target. 


The program is written in FORTRAN G for the IBM 360 computer. 
ITI. APPLICATIONS 
The following applications of this program have been utilized: 


a. Projectile response and trajectory ~ the state of the 
projectile as a function of time for a given situation may be obtained, 


b. Guidance accuracy ~- a Sample distribution of impact locations 
about the. target may be generated, 


c, Footprint analysis ~ the locus of feasible target locations 
for a given combination of zone, QE, enablement time, and environmental 
variables may be obtained, 


d. Parametric studies ~ the effect of varying system, environmental, 
or scenario parameters upon any of the above may be studied. Likewise, 
different systems may be compared under similar conditions. 


IV. SCENARIO DESCRIPTION 


The basic scenario consists of a flat ground plane upon which a 
target is located, a projectile which flies from its enablement point 
onward, presumably to acquire and home into the target, a laser designator, 
and an environment, 


A, Initial Conditions, The simulation begins with the projectile 
in an "acquisition-enabled" state. It has not yet acquired the target but 
is capable of doing so as soon as a sufficiently strong signature is 
received, The positions and velocities of both projectile and target 
are specified; the attitude of the projectile is also specified, 


B. Environment. The target is located upon a flat ground plane 
having a specified altitude with respect to sea level, Other variables 
characterizing ambient conditions include a cloud ceiling, visibility 
range, and wind vector, 
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C. Sequence of Events. The normal sequence of processes modeled is: 
search, acquisition, null, guidance, and impact. Loss of acquisition may 
occur at any time if weak or missing pulses occur or if the target leaves 
the field of view, Reacquisition may commence upon subsequent reappearance 
of the signature, Also, the projectile may be programmed to transition 
from the laser-destgnated mode to a passive infrared-tracking mode, 


D. Noise Sources. The Cypes of noise modeled are spot motion and 
pulse dropout, pitch-yaw gimbal cross-coupling, gyro drift, control surface 
misalignment, and IR detector jitter, The principal contributor to dis- 
persion of impacts in the semi-active laser (SAL) guidance mode is spot 
motion. Gyro drift can be an important source of unguided delivery error 
for gliding projectiles. It is a source of a false steering signal during 
guidance for gliding or non~gliding projectiles, Control surface misalign- 
ment is modeled as another source of a false steering signal during guidance. 
The IR jitter is a white (uncorrelated) angular error internal to the 
detector used in the bifunctional passive IR (BF) mode and is the principal 
contributor to dispersion of impacts in the BF mode of guidance. 


V. FLIGHT EQUATIONS 


A. Inertial, Displacement of the projectile center of mass is CC 
governed by the vector equation 


->> 
where F is the vector sum of aerodynamic, gravitational, and rocket motor 


forces. The optional rocket booster motor is characterized by a rectangular 
thrust-time curve; that is, the thrust is either the nominal specified 
value or it is zero, depending upon the time, Likewise, the fuel mass 
decreases at a constant rate during burn, 
Aerodynamic pitching and yawing of the projectile are approximated 
by a second-order differential equation. 
ate) +2707 Aata=K6 
a aa a 
where 


O& is the angle of attack’ (in yaw or pitch), 


6 is the corresponding control surface deflection, 


6 is the attitude of the velocity vector, 


tT. is a "time constant" varying with the projectile's aerodynamic 
state and depending on its design, / 


K. is a gain also dependent on projectile state and design, 


Coy is the natural pitch-rate damping coefficient of the projectile, 
a constant. 


The implementation of these equations is by simple rectangular 
integration using a time step significantly smaller than the physical 
time constant involved, T_, which is normally of the order of a tenth of 
a second. (The time scales of translational accelerations and velocity- 
orientational accelerations are normally much larger,) The trajectories 
obtained using this simple scheme have compared favorably with those of 
other methods, including 4th-order Runge-Kutta. An earlier version of 
this program was successfully used to reproduce an actual guided trajectory 
(MGP-7) within the precision of the test data. 


The computational scheme for the above equations may be described 
as follows: If we let the subscript j denote the value of a variable 
at time step j, and k denote the value at the succeeding time step k = jt+l, 
then 


a = a h ‘xj 9 
= 
xk xj 0 ad : 
where 
as - Fj/m 


for translational component x, and 


Oo, = O +h Oa 

® _ € 2 e 

a = O, oo h[K,9,/T,, - (a, /T, + 2 Gy%4)/T I] 
~ (6, ~ 64) +h F/T, 


for yaw or pitch angle of attack, The last term above represents an angular 
acceleration due to the malthrust moment of a booster rocket motor. The 
constant h is the integration time step. 


The sign conventions used in these equations are: positive pitch 
angle of attack is nose upward; positive yaw angle of attack is nose left; 


positive control surface deflection is in the sense required to produce a 
positive angle of attack, 


The projectile's roll attitude is fixed. For simplicity, the roll 
attitude is such that pitch and yaw are vertical and horizontal, respectively. 


B, Aerodynamic, The aerodynamic forces are obtained using aero 
coefficient tables, Input to the tables are the mach number, angles of 
attack, and control deflections, The tables are interpolated linearly, 


The normal force coefficients are computed as 


= T\ 6 3 
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where a' is the net angle of attack and M is the mach number, 


The axial force coefficient is computed as one of the following, 
depending on the format of the tables provided by the developer: 
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The speed of sound a and air density € are obtained from a subroutine 
which uses piecewise curve-fits with altitude as the independent variable. 
Then the projectile's mach number M is related to the air speed v as M=v/a. 


The body forces are first computed as 


Fay = Cyy * 4d >» 
Fup = Cup ° Aq , 
Eg 6, ha 


A is the projectile cross section, and v is the air speed. These forces 
are then resolved by rotation into the x,y,z frame of reference, after 
which the gravitational and motor components are added. 


The tables also give K_(M,6) for yaw and pitch and the static margin 
in calibers x » which are used to compute Lm for yaw and pitch as 


8T 
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where I, is the transverse moment of inertia and D is the reference 
diameter of the projectile, 


All of the computations discussed in this section are repeated every 
time step. 


VI. GUIDANCE ANB CONTROL 


A, Acquisition. Two acquisition models are incorporated, a 
deterministic acquisition range model and an energy-flux density model. 


The first is used to determine when acquisition first occurs when 
either (1) the internal mode or (2) the mixed (external position, internal 
intensity) mode of spot motion is selected. The second model is used 
after entering within acquisition range in either of these modes. It is 
also used exclusively for (3) the external spot motion mode, 


The acquisition range model is an equation giving the range at which 
the flux density equals the acquisition threshold for specific values of 
a set of variables, One of these variables, the angle between the line of 
sight and the normal to the designated surface, varies continually with 
each appearance of a spot; thus, the acquisition range is recomputed for 
every spot until it exceeds the actual slant range, 


The flux density model is a more general model relating the flux 
density at the seeker to another set of factors, One of these, the 
scattering cross section, is either a constant defined at acquisition in; 
spot motion modes (1) and (2), or a pulse-to-pulse variable obtained from 
the input cards in mode (3). The flux-density model requires the flux 
density to be both greater than an absolute threshold and not less than a 
specified fraction of the previous flux density to be visible. 


Acquisition is established upon the reception of a given number of 
consecutive pulses. Loss of acquisition occurs if a given number of 
consecutive pulses go unseen, Figure 1 shows the logical sequence of 
acquisition and loss of acquisition. 


B. Detector Characteristics. The detector thresholds, absolute and 
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a. 


relative, have already been mentioned. For a pulse to be received, it 
must meet those criteria and fall within the seeker's instantaneous field 
of view (FOV), which is a circular field, The position and intensity of a 
received pulse determine, by means of a table of transfer characteristics, 
the commanded gyro slew rate and the estimated line-of-sight (LOS) rate, 
Specifically, the detector characteristic is conceived in this program 

as an input-output relationship between two continuous-variable inputs and 
two outputs, The inputs are: (1) the angular difference in degrees in 
pitch (or yaw) between the LOS to the laser spot and the center of the FOV, 
and (2) the signal level of the received pulse in db above the threshold 
level, The outputs are: (1) the estimated LOS rate in degrees per second, 
which is input to the autopilot, and (2) the commanded gyro torquing rate 
in degrees per seond, which is input to the gyro torquing circuit. Moreover, 
the detector characteristic may depend upon two logical conditions: (1) 
whether the detector/autopilot is in the nulling or guiding mode, and (2) 
whether the detector is in the SAL or BF mode, 


These transfer characteristics are entered as input data and are 
interpolated linearly with respect to deviation (apparent angular position 
of the spot on the seeker face) and signal level (db of energy flux density 
above threshold, 10 db = 1 decade), The interpolation scheme permits the 
user to "run out of the table" or extrapolate beyond the ranges of the tabu- 
lated independent variables--it simply assumes the value appropriate to the 
nearest tabulated value of the exceeded variable, (This "plateau" extra- 
polation scheme, incidentally, is also used in the aerodynamic tables, ) 


These tables are assumed to be (1) antisymmetric with respect to null 
(center of FOV) and (2) identical in yaw and pitch channels, Figure 2 
shows a set of hypothetical transfer characteristics, and Table 1 shows the 
form in which these curves would be entered into the program, (This 
example does not represent any actual design,) The arrays are variable- 
dimensioned according to array sizes input. In addition, a second slew 
rate array is provided for the case that torquing behavior is different in 
the gyro-nulling phase (the period between acquisition and commencement of 
guidance) than in the guidance phase, If the BF detector has different 
characteristics than the SAL detector, another complete set of arrays (using 
the common set of deviations) is provided. 


All abscissas for which there is a break point in any of the curves 
(Fig, 2) must be entered, along with the corresponding ordinates for every 
curve. Note that steps in ordinate are permitted by entering the abscissa 
twice, with the ordinates in the order of occurrence. Also note that the 
deviation and signal level 0 must be entered explicitly, 


To the gyro slew rate thus defined by deviation and signal level are 
added the torquer-loop g-bias slew rate and the gyro drift rate. Torquer 
bias acts only during guidance. Drift acts when the gyro is in either the 
track or free mode; in the electrical-cage mode its only effect is to bias 
the gyro slightly off the caged position. 


C. Autopilot. The ZOT.14 autopilot can be used to simulate all the 
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TABLE 1. HYPOTHETICAL DETECTOR CHARACTERISTIC 
(TABULAR FORM) 


DEVIATION (Deg) 


ESTIMATED LOS RATE 


3 3 3 


7 7 5 
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non-gliding systems thus far encountered and the Martin ED gliding systen, 
which keeps the body axis aligned with the "free" gyro axis. The pitch 
channel of the basic autopilot in the guidance phase is represented by 
Figure 3. 


The input \_ in Fig. 3, the estimated LOS rate, is the result of the 
seeker transfer‘ outputs cross-coupled in pitch and yaw as follows: 


p(max)| 7 


where the primes indicate seeker transfer outputs, I is the pitch gimbal 
angle, and r is the cross-coupling coefficient, P 


The input l ts obtained by differencing I at successive time steps, 
and I’ by differencing body-axis and gyro~axis orientation angles. 


The first-order differential equation in the proportional-navigation 
channel (the upper channel in Fig. 3) is rectangularly integrated, as were 
the flight equations: 


6 ) +h (Koh “- 5p - 6 


PNk  °PN4 puy?/™s » 
where 6,,. is the component of deflection due to the proportional navigation 
channel of the autopilot, and 5, is the g-bias deflection, 


The other differential equation, however, typically has a time constant 
T, of about the same order as the time step used, Rather than decrease the 
gine step by an order of magnitude or more, we have assumed that the input 
is constant (with the value calculated above) during the present time step 
and solve the differential equation exactly: 


~h/T 


6 ~ PIKy (1 -e [e7b/T 9}, 


SDk QW] + S¢p 


j 


where 6 is the component of deflection due to the synthetic pitch-rate 
damping channel, and the expressions in brackets are precomputed constants. 


Note in Figure 3 that the resulting deflections are limited in two 
places, 


Different projectiles have not only different parameter values, but also 
different sequences of events in determining the bias 6. and in freeing, 
caging, and tracking the gyro, These differences are dteckheed in Section 


VIT. TARGET MOTION 


A, Target Model, The target is modeled as two distinct geometric 
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points: 


(1) The "spot", which has a position determined jointly by a 
deterministic target maneuver model and a stochastic spot motion model; 
and 


(2) The “target'" or intended impact point, whose position is 
determined solely by the target maneuver model. 


Note that while target (2) is the reference point for miss distances, 
target (1) is what the projectile homes on in the laser—guided mode. 

(There is also a bifunctional mode, in which the projectile may switch from 
tracking target (1) to target (2).) 


B. Maneuver Model. This model provides the capability of four types 
of motion: (1) uniform motion; (2) uniform acceleration parallel to heading; 
(3) uniform turns (constant speed and radius); and (4) sinusoidal zigzags. 
Maneuvers may be programmed in any desired sequence, with the limitation 
that the total number of maneuvers of types 2, 3, and 4 for both targets 
together must not exceed ten, The two targets are programmed independently, 


C, Spot Motion Model, The random motion of the spot about the nominal 
designation point is simulated by either of two means; in either case, it is 
necessary to do some preprocessing, The ERIM target reflectivity model, 
which utilizes a digital-Butterworth-filtered white noise source input to 
obtain a record of designator angular errors, a faceted target model, and 
an assumed designator-target~-seeker geometry, produces a record of apparent 
spot positions in 3-dimensional space. This record may be used to drive 
ZOT.14 directly (external spot motion modes 2 and 3) or its autospectrum 
and spatial covariance matrix may be transmitted to ZOT.14, which creates 
a time series of spot positions having those properties (internal spot 
motion mode 1). 


The internal mode operates in three logical stages: (1) generation 
of three independent gaussian white-noise sources by sampling standard 
normal deviates; (2) recreating the autospectrum by passing the white noise 
components through a digital Butterworth filter; and (3) recreating the 
covariance matrix by performing a linear transformation upon the independent 
filtered noise components. This mode has the advantage of producing an 
indefinitely long record using a small amount of input, 


D, Pulse Dropout. Missing pulses may occur in the real world from 
several causes, such as drifting to a point hidden from the seeker's 
point of view (masking) or out of the seeker's field of view or to an area 
of very low reflectivity, When a missing pulse occurs in the ERIM model, 
zero~scattering cross section and zero position are output for that pulse; 
the ZOT,14 external spot motion mode interprets the zero cross section 
correctly as a missing pulse and the seeker model responds appropriately 
according to the seeker design specifications, In the internal mode, 
ZOT.14 provides an alternating stochastic renewal process of "see" and 
"dropped" time intervals, 
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E. Random Number Blocking. It may be desirable to reproduce identical 
spot position, pulse dropout, and IR jitter sequences for corresponding 
replications of several experiments, This requires that the same initial 
seed values be used in each experiment and that corresponding replications 
"use" the same number of laser and IR spots and pulse dropouts, This second 
requirement is met by requiring all replications to "use" a fixed "block" of 
random numbers regardless of the amount actually required within any one ~ 
replication, : 


Block lengths for spot motion, pulse dropout, and IR jitter are named 
LENGTH, LENGT2, and LENGT3, respectively. These are input parameters and 
must be either all zero or all positive. If positive, they are the allocated 
block sizes; if these lengths are exceeded in any replication, the program 
prints a notification and aborts, If zero, they are internally defined at 
the successful completion of the first replication, allowing for two 
seconds of spot motion and IR jitter more than were actually used in that 
replication and ten cycles of pulse dropout more than were used, In 
practice, the user would make single-replication reference runs for each 
condition of an experiment expected to affect block lengths. The user would 
permit the program to define block lengths for each of these runs and would 
then use the maximum of each length in his experiment, 


This routine operates in all spot motion modes, However, pulse dropout 
blocking operates only in internal mode (1) even though LENGT2 is defined 
for all modes, | 


VIII. FLIGHT TERMINATION é> 


A. Decision Criteria, Termination of a replication occurs upon any 
of the following events: 


(1) Excessive simulated time elapses., The time limit is 1.5 
times initial slant range divided by initial projectile speed, The run is 
aborted. 


(2) A random number block length is exceeded, The run is aborted. 


(3) The projectile achieves closest approach to target (2). 
This point is determined by checking the sign of the dot product of the 
displacement vector from target to projectile with the projectile velocity 
vector. If the product is negative, the projectile is approaching the 
target; if positive, it is receding. The run is aborted if, at closest 
approach, the seeker has never acquired the target; otherwise, the run 
continues with impact scoring and the next replication. 


(4) The projectile impacts the ground. Program action as in 
(3) above. 


B. Scoring. When a replication terminates by criterion (3) or (4) 
and acquisition had occurred, the impact is scored; that is, the position 
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of the point of closest approach is computed relative to target (2), This 


position is given both in the x, y, z simulation coordinates, and in the 
© "impact plane," 
normal to the relative velocity vector, and having yaw (horizontal) and 
pitch (vertical) coordinates, 


The computations proceed as follows: let 


a be the relative location of the projectile at the time 
step preceding the impact decision, 
Vii» Vi» v__ be the relative velocity at that time, 
xF ye 2r 


Yh be the horizontal resultant of that velocity, 
= be the magnitude of that velocity, 


ee oe be the relative point of closest approach, 


r be the miss distance, 


y be the yaw component of miss distance, and 


ey p be the pitch component, 
Then D D 
%, Ve 4 Ver pe Yer V yey a ae 
xX = 
ca 2 
Vv 
Pp 


and similarly for Yon? as by cyclic permutation of coordinates. Then, — 


r= \fx + TZ 
ca Yea ca \ 
Vv ~ Vv 
xr ca zr ca 
y= . » and 
h 
Z 
“VV ry,“ Fy ~VeYv 
: su» ee. Me Ca h_ “ca Go 2r Cay « 
P Vv, Vv 
h p 


given v, and i, + 0. 


IX. REINITIALIZATION 


EF After scoring, the next replication begins with reinitialization of 


a 


a hypothetical plane having origin at target (2), oriented 


the simulation. The unused portions of the previous replication's random 
number blocks are generated and discarded. Then, the projectile and target 
Status are recalled for the reinitialization point, which is one of two 
points: (1) if the target location does not vary between replications and 
gyro drift is deterministic, the reinitialization point is the point at 
which acquisition was first established; otherwise (2) the reinitialization point 
is the original initialization point. When reinitialization is at point (1) 
above, sufficient spot motion is run to align the spot motion block with 

the simulation time of reinitialization, The pulse dropout block is not 
advanced, because the pulse dropout routine functions only after acquisition 
has first been established. 


X. PROJECTILE MODELS 


The various projectile designs accommodated by this program differ in 
the following qualitative respects: 


(1) trajectory -- ballistic or gliding, 


(2) gravity bias -- none, fixed, or adaptive, Bias may be applied 
through the proportional nagivation circuit alone or through the gyro 
torquing circuit as well. 


(3) guidance enablement strategy -- this is the behavior of the 
projectile between acquisition and the commencement of proportional 
navigation guidance, 


(4) rocket booster -- present or absent, ® 
(5) bifunctional tracking -- present or absent, 
(6) behavior upon loss of acquisition, 

Each of these points is discussed in the following paragraphs, 


A, Trajectory. The ballistic trajectory is typical of all but the 
Martin ED CLGP in glide mode. For the TI and Navy projectiles, the gyro 
is electrically caged to the projectile form axis between spinup and 
acquisition, while in the Martin AD case a gyro cant (or lookdown angle) 
was prescribed, The Martin ED ballistic mode is uncanted, like TI and Navy, 
The small lag in the electrically caged gyro may be accounted tor, Lf 
desired, by means of the cant option. The cant is a pitchwise bias applied 
to gyro position relative to the form axis only during the preacquisition 
phase of flight. During ballistic flight, the control surfaces of the 
Navy projectile are in the trail position; all of the Army models have a 
synthetic pitch/yaw rate damping circuit which is active after spinup (thus, 
it is active after acquisition-enable), so that some control surface 
deflection may occur during "ballistic" flight. : 


The Martin ED glide mode projectile is initialized in glide trim; the 
gyro attitude and the deflection and trim angles required for zero transverse 
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acceleration at initialization, given projectile position and velocity, are 
computed by an iterative procedure. The glide mode is an attitude-hold 

mode which "cages" the form axis on a free gyro, in effect the reverse of 

the ballistic situation, ZOT.14 maintains this attitude in an ad hoc fashion 
as follows: at each time step, the bias deflection required to align the 
axis is computed as 


Ss, * = TM/K, 


where 6_ is the required bias, ao is the present angle of attack, I' is the 
present gimbal angle, and K, is the ratio of trim angle of attack to control 
surface deflection. The variable 5, passes through the proportional 
navigation (PN) channel before being applied to the control surfaces, The 
following figure and equation explain the rationale of the above equation: 


gyro axis = desired form axis 


present form axis 


velocity 


a-f 


Cdesired * Ky Sdesired 


B, Gravity Bias, G-bias is present only for the various Martin 
projectiles, For the AD version, it is a fixed value applied through the 
proportional navigation (PN) channel only; for the ED version, it is 
an internally-computed value applied through both PN and gyro torquer 
channels, However, since the computation of the bias is accomplished 
prior to the simulation's initial point, it is entered as a constant in 
the ballistic mode. Any value input in glide mode is ignored during the 
preacquisition glide, since the program computes the initial trim internally, 
The value of bias: in- glide mode is updated as long as the preacquisition 
glide continues, Thereafter, it is fixed at either the last computed 
value or the input value, as selected (there is no actual projectile which 
switches to a preset value -- the Martin ED holds the last computed value). 


Note that G-bias is conceived as a deflection and not as an LOS rate. 
However, during active guidance, the torquer bias channel causes the gyro 
to drift at a fixed rate to introduce an artificial upward LOS rate into 
the PN channel, This is used in the Martin ED projectile to achieve the 
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desired overbias. Since the optimum level of bias during guidance is 
somewhat greater than the one gee which glide requires, this is one way 

to increase the bias at the appropriate time; in fact, this is the way 
Martin has chosen. Also, if acquisition is lost, the bias returns to the 
one-gee level; for all biased models, the deflection bias remains on during 
a lost-acquisition period. 


C. Guidance Enablement. Several strategies are implemented: the LOS 
rate input to the PN channel is shut off until either (1) the gyro has 
first been torqued according to its characteristic to within a nominal 
angular distance from null and a Subsequent delay has then elapsed; or (2) 
the gyro has torqued according to its characteristics for a fixed time. 
The procedure of (1) is accomplished at reacquisition as well; strategy 
(2) permits immediate guidance upon reacquisition regardless of initial LOS 
error. During the nulling period, the projectile may, if specified, be 
steered in a "pursuit" mode in which the control deflection is equal to 
the deflection applied before acquisition (i.e., the glide bias if gliding 
and zero if ballistic) plus a term proportional to the gimbal angle. 

Note that during the nulling period, the gimbal angle rate gives rise to 
a small perturbing control deflection via the Synthetic damping circuit, 
which in effect attributes all gimbal rates to body rates, not gyro rates. 


The purpose of any nulling Strategy (guidance enablement strategy) 
is to allow the gyro-LOS heading error initially existing upon acquisition 
to be eliminated before attempting PN guidance. Without some such 
strategy, the projectile will act upon erroneous steering information. 


D. Rocket Booster. The Navy 5-inch guided projectile has a rocket 
booster, for which two effects are modeled: thrust and torque due to 
misalignment of the thrust vector and center of mass (malthrust). Both of 
these are modeled deterministically. The malthrust moment vector can be 
made to rotate about the form axis at a fixed rate to account for the 
effect of the rolling of the projectile's airframe upon the direction 
of the malthrust. 


E. Bifunctional Tracking. For the study of a projectile which can 
initially acquire a laser spot and Subsequently switch to a passive IR 
Signature, the bifunctional option is provided. The transition occurs 
Whenever (1) slant range to the actual target is less than a deterministic 
transition range and (2) the actual target is within the IR FOV. In this 
model no spot motion (a property of the designator, target, and projectile, 
jointly) is present; however, an internal noise source, IR detector jitter, 
is present. This noise is modeled as a normally-distributed, uncorrelated 
angular error in the detected LOS angles in yaw and pitch. 


F. Behavior upon Loss of Acquisition. Ballistic models recage their 


gyros according to their torquing characteristics; gliding models let 
their gyros go free. Whatever deflection bias was present in guidance 
remains applied. During the decorrelation sequence immediately preceding 
loss of acquisition, gyros are freed and estimated LOS rates hold their 
previous values. 
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XI. PREPROCESSING 


Preprocessing consists of three phases: scenario selection, zoning 
solution, and signature development, 


A, Scenario Selection. In this first phase, the user chooses gun, 
designator, and target locations upon a ground plane, The nature of the 
gun, designator, target and projectile are all specified as well. Finally, 
the environment is characterized by such factors as ground plane altitude, 
visibility, and cloud height, 


B. Zoning Solution. The second phase utilizes the exterior ballistic 
simulation program EXBAL to compute a family of trajectories from which 
a suitable solution for the problem can be taken. 


C. Signature Development, The ERIM target reflectivity simulation 
model generates a record of apparent spot positions as seen from the 
seeker's point of view, 


Input data include: the designator position and direction of approach 
of the projectile, in the target frame of reference; a geometric description 
of the surface of the target vehicle; and the dynamic characteristics of 
designator spot motion. 


The record of spot positions is further processed to compute the 
spectral density and the covariance square root matrix, a matrix which is 
the matrix square root of the covariance matrix, This square root matrix 
is used to transform uncorrelated noise records into properly correlated 
noise records in ZOT,14, Alternatively, the ERIM spot position record may 
be used in ZOT,.14, directly. 


XII, POSTPROCESSING 


The postprocessing program IMPAC may be used for further scoring of 
impacts. This program uses a simplified geometric description of the 
target vehicle and the record of impacts punched onto cards by ZOT.14 to 
score each impact as a hit or a miss, Hits may be further scored by 


computing the impact obliquity (angle between the projectile form axis 
and the normal to the target surface element struck), 


Next page is blank. 
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APPENDIX A 


USER'S GUIDE 


Next page is blank. 
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INPUT 


Input to ZOT.14 consists of a set of at least 15 punched cards. 
Some cards are optional, others required. The card set is assembled 
as follows: 


CARD 1, required, format (8110), contents by field: 


(1) 


(2) 


(3) 
(4) 


(5) 


(6) 


(7) 


(8) 


CARD 


(1) 
(2) 
(3) 


(4) 


MRUN 


MODEL 


NREPS 


NSAMPL 


NPRINT 


INTRVL 


NPUNCH 


MEASUR 


2, required, 
ISEED 
ISEED2 


LENGTH 


LENGT 2 


positive integer -- user's run identification 
number 


integer identifier of aero coefficient data -- 
1 - Martin Marietta AD 

2 - Texas Instruments AD 

3 - Navy 5"/155mm sleeved round, AD 

4 - Martin Marietta ED 


maximum number of replications to be run 


number of integration time steps between designator 
pulses 


print control code 

0 - print final results of replications only 

1, 2 - print also intermediate results at intervals 
as speciffed by INTRVL and print notices of key 
events 


print interval, number of integration time steps 
between printouts 


punch control code 

O - no punch 

1 - punch results of each replication onto a pair of 
cards for IMPAC program 


units of measure control code 

1 - input interpreted as being in metric units (mks) 
2 - input in English units 

format (7110,6X,4I1), contents by field: 

positive odd seed for spot motion generation 


positive odd seed for pulse dropout generation 


block length for spot motion; if zero, LENGTH will 
be defined at end of first replication 


block length for pulse dropout; if LENGTH (above) 


is zero, LENGT2 will be defined at end of first 
replication 
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(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


CARD 3, required, 


(1) 
(2) 
(3) 


(4) 


(5) 
(6) 
(7) 
(8) 


MODESM 


NSL 


MODDEP 


NMODES 


XP 


YP 


ZP 


VZP 


H 


GROUND 


spot motion mode code 

1 - internal generation mode 

2 - mixed external position and internal intensity 
3 - external generation mode 


number of consecutive seen pulses required for 
correlation (acquisition) 


number of consecutive unseen pulses required for 
decorrelation (loss of acquisition) 


number of LOS deviations tabulated in the detector 
characteristic tables (see also CARD 10 and para. VI, 
B) 


number of signal levels in detector characteristic 
tables 


gyro torquing mode dependence code 
1 - gyro torquing characteristic not mode dependent 
2 - torquing characteristic different for nulling and 
guiding phases 
number of distinct detector characteristic sets 
1 - SAL only, or IR detector characteristics are 
identical to SAL 
2 - IR detector characteristics are different from SAL 
format (8F10.0), contents by field: 
“ 
x-coordinate of initial projectile position [m;ft] 


y~-coordinate (altitude above ground) 


z-coordinate (Note: (x,y,z) constitutes a right- 
hand coordinate system 


x-component of initial projectile velocity with 
respect to earth [m/s;ft/s] 


y-component 
z-component 
integration time step [sec] 


ground plane altitude above mean sea level [m;ft] 


” 
The first unit in brackets applies when MEASUR=1, the second when MEASUR=2., 
If only one unit appears, it applies in either case. 
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CARD 4, required, format (8F10.0), contents by field: 


(1) 


(2) 


(3) 
(4) 
(5) 
(6) 


(7) 
(8) 


CARD 5, required; 


(1) 


(2) 
(3) 


(4) 


(5) 


(6) 


Ce, 


(8) 


M 


IB 


KDEL1 
KDEL2 
TDEL 


ZALF 


BIAS2 


CANT /KLDOT 


- FOV(1) 


FOV (2) 


COMLIM 


CONLIM 


KQ 


TAUQ 


DLNULL 


WAIT 


projectile mass, exclusive of fuel if rocket- 

assisted [kg;lbm] (Note: all variables listed 
under F-formats are real-valued, regardless of 
initial letter of name.) 


projectile transverse moment of inertia [kgem’; 
slug*ft2] 


autopilot gain for SAL guidance [sec] Ks 
il 

autopilot gain for IR guidance [sec] Ks 

2 


autopilot time constant Ts [sec] 


projectile natural pitch rate damping coefficient 
[dimensionless] 


constant gravity-compensation deflection bias [deg] 
for a ballistic-mode projectile, CANT is 

the gyro search cant angle, the angle downward from 
the projectile axis to the gyro axis [deg]; for a 
glide-mode projectile, KLDOT is the torquer-loop 
bias gain Ks [sec71] 

Format (8F10.0); contents by field: 


instantaneous field of view for SAL, center of FOV 
to edge [deg] 


field of view for IR 


limit on proportional navigation deflection commands 
[deg] 


limit on control surface deflection [deg] 
synthetic pitch-rate damping gain [sec] 


synthetic damping time constant -- must be positéve 
[sec] 


null zone -- line-of-sight deviation must be less than 
DLNULL to accomplish null 


post null delay -- time during which null continues 
after achieving DLNULL [sec] 


2/ 


CARD 6, required, format (8F10.0), contents by field: 


(1) 
(2) 
(3) 
(4) 
(5) 


(6) 
(7) 


(8) 


CARD 7, required, 


(1) 


(2) 


(3) 


FUELMS 


SIMPLS 


THRUST 


TBSTRT 


YLEVER 


PLEVER 


GDMAG 


GDDIR 


CSMA1 


CSMA2 


ROLRAT 


mass of fuel of rocket motor [kg; lbm] 

specific impulse of rocket [N*sec/kg; 1bf*sec/1bm] 
rocket thrust [N;1bf] 

Clock time (from launch) at burn start [sec] 


rocket malthrust lever arm component creating 
positive yaw moment [m;in] 


malthrust lever arm, pitch 


if positive, deterministic gyro drift rate; if negative, 
absolute value is max random gyro drift rate [deg/sec] 


if GDMAG pesitive, GDDIR is gyro drift direction CCW 
from the right [deg]; if GDMAG negative, GDDIR is 
seed number for gyro drift (positive odd) 


format (8F10.0), contents by field: 


if CSMA2 < 100.0, CSMA1 is the misalignment of the 
pitch control surface; if CSMA2 > 100.0, CSMA1 is 
the standard deviation of random control surface 
misalignment [deg] 


if < 100.0, control surface misalignment is 
deterministic and CSMA2 is the misalignment of the 
yaw control surface [deg]; if > 100.0, misalignment 
is random and CSMA2 is the seed number (positive 
odd) 


angular rate at which the malthrust moment rotates 
about the form axis [deg/sec] 


Remaining fields are not used, 


CARD 8, required, format (8F10.0), contents by field: 


(1) 
(2) 


(3) 


(4) 


GAMAX 


XCUPLM 


TIMEO 


AVGO 


max permissible gimbal angle [deg] 


gimbal cross-coupling coefficient at max gimbal 
angle [dimensionless] 


clock time (from launch) at initialization of 
simulation [sec] 


mean length of a dropped-pulse string [sec] 
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(5) AVG1 mean length of a visible-pulse string [sec] 
(6)  VWIND wind speed [m/s;ft/s] 


(7) THWIND azimuthal direction from which wind is blowing, 
measured CCW from +z axis [deg] 


(8) GLIDON control code number, 3 digits followed by decimal: 
ist digit = 0 for ballistic, 1 for glide mode 
2nd digit = 0 for g-bias = BIAS2, 1 for g-bias = 
computed glide bias, during the guided phase of 
flight 
3rd digit = 0 for null required only once, = 1 
for null required on all acquisitions and 
reacquisitions 


CARD 9, required, format (8F10.0), contents by field: 


(1) 12 at element of spot motion transform matrix [m;ft] 

(2) P12 

(3) P13 

(4) P22 

(5) P23 

(6) P33 

(7) FREQ corner frequency of spot motion power spectrum [hz] 
(8) DR designation range [km;mi] 


CARD 10, at least 1 card required, format (20F4,.0), contents: first 
NID* fields - vector of tabulated LOS deviations [deg], next NSL fields - 
vector of tabulated signal levels [db], next NTDxNSLxMODDEPxNMODES fields 
- table of estimated LOS rates as follows: 
((((TMD(I,J,K,L) ,I=1,NTD) ,J=1,NSL) ,K=1,MODDEP), L=1,NMODES) 
[deg/sec] 
next NIDXNSLXNMODES fields - table of commanded gyro torquing rates as 
follows: 
(((TSR(I,J,L) ,I=1,NTD) ,J=1,NSL) ,L=1,NMODES) [deg/sec] 


CARD 11, required, format (8F10.0), contents by field: 


(1)  THVCN angle between outward normal to designated target 
surface and the vertical [deg] 

(2)  PHVCN azimuthal angle of outward normal to designated target 
surface [deg]. Angle measured counterclockwise fron 
Z- axis 


* 
NID, NSL, MODDEP, and NMODES are defined in CARD 2. 


Zo 


(3) LP laser pulse energy [J] 


(4) THR detector threshold in units of [alos ent ies 
(5) REFL target surface reflectivity [dimensionless] 
(6) VR visibility range [km;mi] | ‘ 
(7) CODE atmospheric water vapor code: 
2. for 7.84 mm/km P 


3. for 1.87 mn/km 
(8)  CEIL cloud ceiling height above ground [m:ft] 


CARD 12, required, format (8F10.0), contents by field: 


(1)  BMDIVG laser beam divergence [mrad] 

(2) DYRANG detector lookdown instantaneous dynamic range [ db ] 

(3) PRANGE visibility range due to particulate scattering 
only [km;mi] 

(4) KAH attitude-hold gain 

(5)  SERRIR standard deviation of white angular noise in passive 
IR detector [mrad] co 

(6) SEEDIR seed for IR jitter (positive odd) 

(7) RLNGT3 block length for IR jitter -- same rules as LENGT2 


(8) not used 
CARD 13, 2 cards required, format (8F10.0) -- 


CARD 13-1, contents by field: 


(1) X(1) x-coordinate of initial reference spot location [msft] 

(2) Y(1) y-coordinate 

(3) Z(1) z-coordinate 

(4)  VEL(1) initial speed of reference spot location [m/s;ft/s] 


(5)  THETA(1) initial azimuthal heading of ref spot location [deg] 


(6) XINCR x-wise increment in location of targets between 
replications [m:ft] 


(7) ZINCR z-wise increment in location 
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(8) 


not used 


CARD 13-2, contents by field: 


(1) 


(2) 
(3) 
(4) 
(5) 
(6) 


(7) 


(8) 


X(2) 


¥(2) 
Z(2) 

VEL (2) 
THETA (2) 


RSWICH 


XINCH 


not used 


x~coordinate of initial location of desired impact 
point [m;ft] 


y~coordinate 

z-coordinate 

initial speed of desired impact point [m/s;ft/s] 
initial azimuthal heading [deg] 


slant range within which projectile may switch to 
passive IR tracking mode [m;ft] 


position units override code -- if not zero, then the 
positions of both targets are interpreted as being 
in inches 


CARD 14, required, format (8110), contents by field: 


(1) 


(2) 


MANVRS (1) 


MANVRS (2) 


number of programmed evasive maneuvers 
for target 1 (the designator signature) 


number of evasive maneuvers for target 2 (the desired 
impact point) 


Remaining fields are not used, 


CARD 15, optional, format (8F10.0) -- 


The number of cards equals the sum MANVRS (1) + MANVRS(2). 
Cards are in the order of execution of maneuver for target 1, followed 
by order of execution for target 2, 


Contents by field: 


(1) duration of target's straight run preceding the maneuver [sec] 


(2) maneuver code: 


2; 


4, 


= acceleration/deceleration 
3, ™ constant speed/radius turn 
= sinusoidal zigzag 


(3) for code 2, absolute value of acceleration Faas | 
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for code 3, angle of turn, positive left (driver's perspective) 
[deg] 
for code 4, number of complete sine waves 


(4) for code 2, final speed [m/s] 
for code 3, turning radius [m] 
for code 4, not used “ 


Remaining fields are not used. 


CARD 16, optional, format (I10,4F10.0) -- 
Used only if MODESM = 2 or 3; then 1 card required,. 


Contents by field: 

(1)  NSPOTS number of spots in external record 

(2) CFACTR multiplicative correction factor to be applied to 
scattering cross-sections of external record, to 


allow for changes in reflectivity 


(3) DESGPT (1) distance forward from vehicle trailing edge to 
mean spot position, times -1 [in] 


(4)  DESGPT(2) height of mean spot position above ground [in] 


(5) DESGPT (3) distance from centerline leftward to mean spot 
location [in] 


Remaining fields are not used. 


CARD 17, optional, format (20X,E10,3,3F10.0) --~ 
Used only if MODESM = 2 or 3; then NSPOTS cards required, 


Contents by field: 
(1) BRITE scattering cross~section ies ieeer 


(2) YV(1) distance forward from vehicle trailing edge to spot 
position [in] 


(3)  YvV(2) distance from centerline leftward (from driver's 
perspective) to spot position [in] 


(4) YV(3) height of spot position [in] 


Card decks thus assembled may be stacked in series; the last deck 
should be followed by a blank card, 
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i, NOTES ON DATA SETUPS FOR VARIOUS APPLICATIONS 


(1) Post-Enablement Phase of Unguided Trajectory, 


Usually, all noise sources are zeroed a To prevent acquisition, 

do not zero laser energy, but rather set ceiling to zero. To prevent 
° premature termination, place the targets considerably farther downrange 
than the projectile is expected to fly, Select one replication, 
NPRINT = 2, and an appropriate print interval (usually the equivalent 
of 0,5 to 2.0 seconds), 


(2) Reference Guided Trajectory, 


Usually, noise sources are zeroed out. Select other parameters as 
appropriate to the scenario, Usually, NREPS=1, NPRINT=2, INTRVL as 
required for a 0.d- to 0.5-sec print interval, 


NOTE: If the time at which the projectile breaks cloud or comes 
within acquisition range is significantly later than the enablement 
time, and if the trajectory dispersion due to random gyro drift is 
negligible or not Present, then much computer time can be saved by 
using a point just prior to breaking cloud or coming within acquisition 
range as the initial point in further runs studying the same scenario, 


(3) Projectile Response, 


{Pp These are usually runs of one to five replications. Some error or 
noise source of interest has been added to the previously noise-free 
scenario of a reference guided trajectory, and the behavior of the 
projectile is to be studied, The print interval ts smaller than 
previously (unless the noise source is of the round-to-round type), 
perhaps equal to or even less than the pulse interval, Any number 
of modeled noise sources may be added in until the complete scenario 
of interest is simulated, 


(4) Guidance Accuracy, 


Guidance accuracy simulations usually account for all pertinent error 
sources, The print level is NPRINT=0, INTRYL=0, and the punch option 
NPUNCH=1 is taken if postprocessing with IMPAC is to be performed, 

The number of replications NREPS is selected with due consideration of 
the precision required in the estimates of guidance accuracy statistics; 
standard statistical techniques give the standard error of an estimate 
as a function of sample size, 


(5) Footprint Analysis, 


Footprints are generated by means of a series of noise-free runs in 


0 — Fi 
Valid seed numbers must be provided for spot motion and pulse dropout, even 
if these error sources are not used, Also, a valid seed number must be 
provided for IR jitter if the passive tracking mode is used. 
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which the target-location incrementing option is utilized. A typical 
starting procedure would be as follows: 
21 replications, incrementing the (stationary) target by 50-m 
steps from the ballistic/pseudoballistic impact point toward 
the gun (uprange)., 


21 replications, incrementing in the same manner from the PBIP 
crossrange., 


21 replications, incrementing downrange, 


Presumably, at some point either the projectile misses the target by 
more than some arbitrary distance (say 5 feet or 1.5m) or the projectile 
is unable to acquire the target, The interval between two successive 
replications (a hit and a miss) plotted on a ground-plane grid, is 

one data "point" defining the boundary of the footprint. 


If a run ends before finding a boundary point, the user sets up 
a new run picking up where the previous one left off, This is commonly 
required in locating the stretch point (farthest downrange point). 


After these three initial points have been located, the user 
can use additional crosstrack~stepping runs to fill in the gaps in 
the boundary to any precision desired, 


(6) Parametric Studies, 


Studies of the effect of varying parameters may take the form of any 
of the above as appropriate. 


OUTPUT 
(1) Printout, 


Printout consists of three parts: data, intermediate results, and 
final results, 


(a) Data section 


The printout begins with approximately three pages repeating 
the input data in a readable format. 


(b) Intermediate results section 


This part of the printout is controlled by NPRINT and INTRVL. 
The intermediate results may consist of only a line indicating 
where the projectile impacted relative to the target, or it may 
also contain notices of events such as acquisition/loss of 
acquisition, pulse dropout/reappearance, target evasive 
maneuver, etc,, and notices of the projectile state at specified 
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time intervals. These results are in order of replication, 
and within a replication, in chronological order, 


(c) Final results section 


This part is printed if multiple replications are specified and 
all replications were completed, It consists of a list of 
guidance accuracy statistics, including ordered miss distances, 


(2) Punched Card Output, 


This output is used as input to the IMPAC program, It consists of 
two cards per replication containing the relative position, the 
velocity, and the attitude of the projectile at closest approach to 
the target. 


Next page is blank 
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PROGRAM LISTING 
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